MPI-Ph/93-40 
Saclay T/93/063 
June 1993 



Finite size scaling analysis 
of intermittency moments 
in the two dimensional Ising model 



Z. Burda and K. Zalewski 

Institute of Physics, Jagellonian University 
Reymonta 4, PL-30 059 Cracow, Poland 

R. Peschanski 

Service de Physique Theorique 
Centre d'Etudes de Saclay, F-91191 Gif-sur-Yvette, Cedex, France 

J. Wosiek 

Max Planck Institut fiir Physik 
W-8000 Munich 40, P.O. Box 401212, Germany 



Abstract 

Finite size scaling is shown to work very well for the block variables used 
in intermittency studies on a 2-d Ising lattice. The intermittency exponents so 
derived exhibit the expected relations to the magnetic critical exponent of the 
model. 
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Introduction. Intermittency studies have developed into an active branch of 
particle physics. In spite of much effort, however, the origin of intermittency in 
multiparticle production processes is still controversial. In the original papers |I| 
two possible mechanisms have been pointed out. Intermittency could be: a result of 
particle production by long random cascades, or a signal of a phase transition. It 
was soon suggested (for early reference cf. 0) that it also could result from suitably 
conspiring ordinary processes like resonance production and decay, Bose-Einstein 
interference etc. Standard QCD cascade was also shown to lead to intermittency 
at very high energies 0. In order to shed some light on possible mechanisms, it is 
instructive to consider solvable systems where intermittency is present and its origin 
is unambiguous. 

It has been pointed out by one of us |Q that intermittency should be observable 
at the phase transition point of the two-dimensional (2d) Ising model. Indeed the 
transition is of second order, hence the fluctuations are selfsimilar at all length scales. 
This is usually the condition for the occurrence of intermittency. In order to define 
it quantitatively, the L x L Ising lattice is divided into cells of size / x I. For each 
cell a block variable ki is defined, and the thermodynamical averages of moments of 
ki are computed using standard Monte Carlo methods. If the moments Mp show a 
power dependence on the size of a cell, i.e. if 

logMp = -Aplog/ + Cp, (1) 

where Xp and Cp do not depend on /, the system is said to be intermittent with the 
intermittency exponents Xp. 

Such a formulation of the problem bears a striking resemblance to the renormal- 
ization group (RG) approach used to study critical phenomena . Indeed, employing 
results of the RG analysis, Satz derived a relation between intermittency expo- 
nents and the known critical magnetic exponent 6 of the Ising model. Let us choose 
as the block variable 

ki = \si\, (2) 

where si = |7 denotes the average spin in an / x / cell Ci. Satz used si 

without the absolute value, which leads to some difficulties, but the idea is the same. 
According to the RG [0 , when calculating thermodynamical averages, the dynamical 
variables si can be replaced by 

~si = I'^si, (3) 

where the new variables si = ±1 just like the original spins. The constant x is known 
as the anomalous dimension of the spin field. For the two-dimensional Ising model x 
can be expressed by the critical magnetic exponent 6 = 15 (cf. P). One finds 

X = = - (4) 

6 + 18 ^ ^ 
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Using Eq.(^) one obtains the power dependence on the block size 

Coefficient g depends in general on all independent ratios of the dimensional pa- 
rameters entering the problem. In particular the dependence on temperature enters 
through the ratio of the infinite volume correlation length to any other relevant 
scale parameter. Assuming that g{0, L/C,{l3)) 7^ 0, the above finite size scaling Ansatz 
gives Eq.(|l]), for 1 << / << L with the intermittency exponents 

K = ^P- (6) 

Let us stress that Eq.(|D represents only the leading, in /, result of the RG analysis. 
In general nonleading terms may affect this simple behaviour. 

A convincing numerical confirmation of Eq.@ was, however, more difficult. In 
Ref. 0] the number of "up" spins in a cell |(1 + s/) was chosen as a block variable. In 
the vicinity of the critical temperature the increase of the normalized moments with 
decreasing cell size was demonstrated qualitatively. Later P, || it was pointed out that 
Z2 symmetric observables may be more convenient for quantitative determination of 
the intermittency indices, since they are insensitive to fluctuations between the two 
ordered phases. Such transitions, which occur in finite systems, may mask the true 
critical fluctuations. 

The authors of Refs [|, H suggested the block variable |(1 + Si), where 

Si = si signisL), (7) 

They also choose to work at the quasicritical temperature where, at given L, the 
magnetic susceptibility is the largest. At ffist they did not get a clear signal for 
the power dependence [||, but an improved analysis uncovered the selfsimilar 
behaviour; however, rather surprisingly, with intermittency indices derived from the 
percolation critical exponent and not from the magnetic one. A possible explanation 
of this result was suggested in [0. Moments of this particular block variables are 
linear combinations of averages of various powers of 5*;. Assuming for each power of 
Si a scaling law analogous to (|l]), one indeed finds, for the lattice sizes considered, 
the effective exponents similar to the percolation ones. Whether this is a numerical 
coincidence, or evidence for a genuine relation to percolation, is an open problem 
[p!T| , p!2| . Very recently Leroyer [|13| used Si as the block variable and for (3 = (3c = 
|log(l + \/2) got a clear con&mation of the power behaviour with intermittency 
exponents derived from the magnetic one. We shall comment on his results later. 

Results. The main points of the present work are: the choice of as the block 
variable, and a consistent finite size scaling (FSS) analysis of our Monte Carlo data. 
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In the FSS approach the problem of choosing the best (3 is less severe, since one is 
studying the dependence on L and /5 in a small scaling window around (3c- In the 2d 
Ising model the critical exponent v is equal to 1, hence the FSS variable is 

/\(3 = (3-(3c. (8) 

Our runs were done for four values of y (0.0,0.105,0.21,0.42), i.e. in the range 
which corresponds to (3 between (3c and the pseudocritical one for given L. We find 
numerically that the results depend weakly on the choice of y in this region. 



Figure 1: Dependence of the second moment, for y = 0.105, on the cell size I. 
Continuous lines correspond to fixed (see text). The dashed line connects points 
for L = 1024. 



MC simulations were performed using the Swendsen-Wang algorithm |14 . The 



values L = 64, 128, 256, 512, 1024 and A <l < L were used. For each lattice size about 
10^ measurements were done. In order to reduce autocorrelations a measurement was 
taken every sixth sweep, for smaller lattices, and every eight sweep for larger ones. We 
took advantage of the very slow growth of the autocorrelation time in the Swendsen- 
Wang method. The critical exponent characterizing this growth is estimated to be 
about .35. For each new lattice the first 10^ sweeps were discarded to allow the 
system to reach equilibrium. Statistical errors of our MC results are always less than 
the sizes of the symbols used to denote points on the figures. At fixed l/L and y, to a 
very good accuracy, the dependences of log < \'si\^ > on log / are found to be linear in 
all the regions considered. The slopes agree with predictions (^) within small errors, 
which can be ascribed to the finite L corrections. 
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In Fig.l we show as an example the dependence of log < > on log / for fixed 
y = 0.105 and fixed l/L = 2~™, m = 0, 1, . . . , 6. The dotted line joins the points with 
fixed L = 1024. The usual approach was to use such lines for I « L to find the 
slopes. Clearly, the slopes of lines evaluated at constant are much better defined 
in accordance with Eq.^). The slopes for the six lines in the figure range from 0.239 
to 0.272 to be compared with j obtained from the magnetic exponent, Eq.(P), and 
II ~ .104 resulting from the percolation exponent. Certainly the magnetic exponent 
is preferred. In Fig. 2 we show the results for y = 0.42, l/L = 1/8 and p = 1,2, 3, 4, 5. 
The lines are normalized at / = 32 and have slopes p/8. The agreement is quite 
satisfactory. 



Figure 2: Dependence of the moments p = 1, . . . , 5 for y = 0.42 on the cell size /. 
The lines are calculated from log2Mp(/) = log2Mp(32) — |(log2/ — 5). 

The normalized moments [|l5l are of little interest for the present choice of the 
block variable, because from Eq.(|^) 

F„ = ^ l^'l^ ^ ~ /°, l/L = const, (9) 

where ~ means equal up to a scale independent factor. At fixed L there could be some 
I dependence from the argument of the scale invariant coefficients. However, since 
our numerical results show that these coefficients tend to nonzero limits as l/L 0, 
there are no relations analogous to Eq.® with non-zero Xp for I « L. 

It is a property of the Ising model, and in fact of many field theories, that 
the first power of a field has a nontrivial anomalous dimension. In the multiparticle 
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language single particle inclusive distributions are not smooth. The anomalous di- 
mension of the n-th power of the field is in our case simply n times the anomalous 
dimension of the field itself. Thus the anomalous dimensions in the numerator and 
denominator of Eq.(^) cancel. Fortunately in order to study selfsimilarity or inter- 
mittency it is enough to normalize moments, Eq.(^, by appropriate powers of the 



volume of a cell which looses no information about critical indices. 

There is also no motivation to use the factorial or binomial [0 moments, 
since, although we always consider the I « L limit, one should keep 1 << / at the 
same time, i.e. the "number of particles" in a cell must be large. This means that 
the bias introduced by the noise is not important, and tricks invented to filter it out 
are not necessary. 

We discuss now moments of the block variable S"/, Eq.(|^. For p even and/or 
I = L we have = \si\^. Thus only the moments with p odd and / < L require a 
separate discussion. As we will see shortly, these moments have a different scaling 
behaviour, which is also confirmed by our data. The main difference with the previous 
block variable is that, due to the translational invariance and additivity, the first 
moment of Si is independent of the cell size [|. Indeed we have 

<Si>=«Si»=<Sl>--L-^^\ (10) 

Where the double brackets denote average over cells and over the ensemble, and ~ 
means, as usual, equal up to a scale invariant coefficient. Assuming that the RG 
result, Eq.(|D, holds for the new variable 

<S^>^l-P/^<Si>, p-odd, (11) 

where we have used Sf = Si for odd p. By definition we have Sl = \sl\ = 1- 
Equations (p!OD and (0) imply 



<Si>^i-] , (12) 



while for the normalized moments, studied also by Leroyer |]T3| one gets 



Fp^^j^, = [^j h{-,y), p-odd, (13) 

where the scale invariant function has been written explicitly. Again, for h{z,y) 
regular and non-vanishing at z = 0, the above formula gives the power behaviour 
with the intermittency exponents 

— p — 1 

Xp = , p — odd. (14) 



^One may say that Si has zero anomalous dimension, however any higher power does not. 
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The behaviour (0) has been numerically confirmed by Ref.[|T3]. On the other hand 
for the unnormalized moments we get 

<sr>=rVL-lif(l y). (15) 

At fixed L » I these moments scale as < > while for / L they coincide with 

< >. Both these features are clearly seen in Fig. 3, where log < > (p = 1, 3, 5) 
and log < > (p = 1, 2, 3, 4, 5) are plotted versus log / at fixed L = 512. Note that 
also < Si > for / << L is parallel to < \si\^ >= 1. 



Figure 3: Dependence of the moments, for y = 0.21, on the cell size /. The dashed 
lines correspond to the block variable \si\ for p = 1,2, 3, 4, 5 and solid lines to Si for 
p = 1,3, 5. Note that the solid lines are parallel to the dashed ones with p ^ p — 1. 

Summary. It is seen that the intermittency of the Ising model is yet another 
consequence of the well known scale invariance of the system in the critical region. 
The renormalization group formula, together with the finite size scaling (^), give a 
good description of this phenomenon in the vicinity of f3c (0 < LA(3 < .42), for 
lattice sizes L = 64 to 1024 and cell sizes / = 4 to L. The power behaviour, Eq.([^), 
was confirmed with the intermittency exponents determined by the magnetic critical 
exponent 5 when \si\ or 5*/ are used as block variables. Moments of both variables 
were shown to obey simple relations which are also confirmed by our data. On the 
other hand, it is still an open question whether it is possible to find a block variable, 
which would yield intermittency exponents related to the percolation exponent. In 



particular the proposal [T^ to consider the largest connected cluster in a cell remains 
to be studied. 
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We would like to stress a few points which considerably simplify both the the- 
oretical and the numerical studies. 

• The absolute value of the resultant spin of a cell \si\ is Z2 symmetric and seems to 
be a convenient choice for the block variable. 

• An advantage of the finite size scaling approach is that instead of selecting a par- 
ticular value of P one can study the simultaneous dependence on P and L by fixing 
y = LAP at several not too large values. This is both simpler and less restrictive. 

• The finite size scaling analysis at fixed l/L is simpler conceptually, and gives better 
defined slopes, than the analysis at fixed L. This implies that (nearly) exact power 
behaviour (|l]) at fixed L can be expected only at / << L (keeping still !<</). 

K. Z. thanks for hospitality the Service de Physique Theorique de Saclay, where 
this work was started, and the Max-Planck-Institute fiir Physik in Miinchen, where 
it has been completed. 
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